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Abstract 

I The statistical properties of heavy particle trajectories in high Reynolds numbers turbulent 

f-H . flows are analyzed. Dimensional analysis assuming Kolmogorov scaling is compared with the 

result of numerical simulation using a synthetic turbulence advecting field. The non-Markovian 
nature of the fluid velocity statistics along the solid particle trajectories is put into evidence, and 
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I. Introduction 

A heavy particle in a turbulent flow moves along a trajectory, which does not coincide with 
that of the fluid parcel where it lied initially 0]. This phenomenon is controlled by the particle 
trajectory relaxation time rs, which, in general, is a function of the difference between the velocities 

V and u respectively of the particle and of the fluid at the particle position. To lowest order in 
the particle Reynolds number Rcp, ts is given by the Stokes time. (For a spherical particle, the 
particle Reynolds number and the Stokes time are defined respectively as Rep = |u — v|a/j/ and 
^■5 ~ ^rl"'^ ^ p/ Pol with p and po the particle and fluid density, v the fluid kinematic viscosity and 
a the particle diameter). Divergence of fluid and particle trajectories receive contribution both 
from external forces (e.g. gravity) acting on the particle and acceleration along a fluid trajectory; 

V obeys therefore, to lowest order in Rcp and po/p, the linear equation: 

v = r-i(u-v)+g (1) 

with (p — po)g the force acting on the particle. Accelerations produced by non-uniformity or non- 
stationarity of the mean flow may be reabsorbed into the term g, after redefinition of u and v, 
through a Galilean shift to the reference frame of the mean flow. 

Starting from the work of Csanady , much of the interest has been in the role of gravity, with 
the conclusion that its main effect on turbulent diffusion is a renormalization of the correlation 
time. A particle, in the presence of an external force, will drift with respect to the fluid with an 
average velocity Tsg, a phenomenon known under the name of trajectory crossing, and this will 
cause a shorter time of permanence of a particle in a correlated region and a Lagrangian correlation 
time shorter than that of a fluid parcel. 

This renormalization has been used for the derivation of Lagrangian models for heavy particles 
dispersion j^], but it is not clear whether a Markovian hypothesis for the velocity is indeed ap- 
plicable to the case of a heavy particle. This problem, associated with the small turbulent scales 
being seen as frozen by the drifting particle, was not present in who dealt essentially with 

a turbulent flow with only large scales present. The same turbulent regime was considered in 
based on an expansion around the case of an infinite inertia particle. In the Eulerian approach of 
, the same problem was avoided working in the limit of Stokes time shorter than the Kolmogorov 
time, which is the same assumption on which the classical theory by Tchen is based 

A different approach is to model fluid and solid particle trajectories on the same ground. The 
result of this line of reasoning is a class of Lagrangian models, of which, one of the prototypes is 



2 



described in |Q. |To|. The basic idea is to couple Eqn. (1) (or its more complex counterpart valid 
for arbitrary Rcp) with the equations for the solid particle trajectory and with the model equation 
for the trajectory of the fluid parcel where the particle lied initially. In high Reynolds numbers, 
one resorts typically to models which are Markovian in the Lagrangian velocity of the fluid parcel 
u^{t) (see e.g. In this regime, the model equations would read, in the absence of a mean 

flow: 



where a(u'',x'') is of the order of /ti^{'x.^), with tl{x.^) the turbulent integral time, is 
the mean turbulent dissipation and Cq is an empirical constant, with values lying between 2 and 
7 depending on the kind of flow x'' is the coordinate of the fluid parcel for which initially: 
x'' = X. After a time A, whose determination is part of the modelling effort, a new fluid parcel is 
considered, thus setting again x'^ = x and the process is repeated; in the meanwhile, when x'' ^ x, 
the problem is to calculate from u'' the fluid velocity u at the solid particle position x. 

This approach is able in principle to deal with both regimes of strong external forcing and 
dominant trajectory crossing, and of random trajectory separation produced solely by finite inertia 
and particles being unable to follow the rapid turbulent fluctuations. Obtaining the fluid velocity 
at the particle position turns out to be the difficult part of the problem. In j9[ 0, the difference 
u — u'' was calculated only at the discrete times t„ — nA, with the values of u — u'' at different 
tn assumed uncorrelated. However, although there have been other attempts in the direction of 
Markov modelling [|l3j, the statistics for u appears to be non-Mar kovian, in general, both in time 
and in the separation x — x''. Only recently p^ , |l5j, has the importance of anomalous scaling 
for the velocity along solid particle trajectories been recognized. However, the mechanism for 
its production was only partially discussed and no prediction for the dependence of the scaling 
exponent on the flow parameters was provided. 

The problems associated with the non-Markovian nature of the fluid velocity u are further 
complicated in the case of weak external forces, with turbulent fluctuations and inertia dominating 
the process of trajectory separation. Analysis of this problem and of its consequence in terms of 
constraints on the possible form of turbulent transport models will be the subject of this paper. 



One can see a turbulent flow as a superposition of vortices of different size r, which, at high 
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II. Scaling analysis 
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Reynolds numbers, obey Kolmogorov scaling. In terms of the space (Eulerian) longitudinal struc- 
ture function, this means, for r ^ L, with L of the order of the size of the largest eddies: 

ul = ([<(x + r,t)- <(x, t)r) = CK{er) % = ul (^) ' (3) 

with u^(x, the Eulerian velocity field, Ck — 2 the Kolmogorov constant and m| = {[uf]'^)- In 
terms of the time (Lagrangian) structure function, for t <^ tl'- 

{[uUt)-uUO)f) = Coet = ul- (4) 

tl 

with Tl of the order of the Lagrangian correlation time T/, : 

poo 

n = uf dt{uUtX{o)). (5) 







The quantities L, tl and ul identify integral length, time and velocity scales of the flow, and, from 
Eqns. (3-4), an eddy turnover time can be defined in the standard way: 



r.=r.(9^ (6) 



Consider a solid particle, which at time t = is in xq moving with speed vq and to fix the ideas 
imagine first that no external forces are present. The trajectory separation will be given by: 



Ax(t) = x(t) - x^(t) = /* dt'[y{t') - u^(t'|xo, 0)] 
Jo 



(7) 



where u''(t|xo, 0) is the fluid velocity of the fluid parcel, which at time t = was in xq. 

Asymptotically, for ts /TLTTr^/rs 0, where rj = viei is the Kolmogorov length and Tt, is 
the associated Kolmogorov time, we can distinguish three phases in the separation process. The 
central role is played by the velocity and length scales of the vortices with life-time rs, namely, for 
Tl > TS-- 

Us = ul{tsItl)^, S = L{ts/tl)^ = usTs; (8) 

Ts and Us will appear to be the time and velocity scale of the velocity difference u — v. 

If uo = u(0) = u^(0|xo,0) is the initial fluid velocity and |vo — uo| <C us, there will be an 
initial phase, in which the particle sees the turbulent fleld as frozen and the separation process is 
ballistic: 

Ax(t) - (vo - Uo)t. (9) 

In fact, if Ax{t) is so small that vortices of size Ax{t) have characteristic velocity WAx(t) ^ |vo— uo|, 
it will also be |^p^*^| <^ T^x{t) i i-e- the crossing time of the vortex by the solid particle will be 



much shorter than the eddy turnover time of that vortex. The contribution from this vortex to 
Aa; will thus be u^x{t) Aa;(f) and then Ax{t) ~ |uo — vo|f. This phase will last until 

WAa;(t) ~ |uo — Vol, when the contribution from the vortices begins to be important. For the 
Kolmogorov spectrum described by Eqn. (2), this leads to a cross-over time: 

TB^C~h-'\vo-Uof, (10) 

and for |vo — Uo| ^ us, it will also be tb <^ts- 

In the next phase, for tb <^t <^ts, the separation process will receive two contributions: one 
from the velocity difference across a vortex of size Ax{t), the other from the fact that the solid 
particle is unable to follow the fluid on time scales much shorter that rg. Prom the relaxation 
equation nature of Eqn. (1), these contributions are both of the same order, and the end result is 
that the particles will separate with a relative velocity of the order of UAx(t)- This leads to a process 
of Richardson diffusion {Ax^{t)) oc t^, although with a coeflacient different from that occurring in 
the case of a pair of fluid particles. Now, for most choices of initial conditions: |vo — uo| ^ us, 
the ballistic phase will end at t ^ tb ''^ Tg and the Richardson diffusion phase will be absent. The 
same will occur when |vo — Uo| > us; the ballistic phase will end when the solid particle trajectory 
relaxes, which happens, again, when tb ~ rg. 

In the end phase, for i » rg, the solid particle will deviate from the fluid parcel it crosses at 
the time t, with a rate much smaller than the rate of separation from the fluid parcel it crossed at 
time zero. In fact, from Eqn. (1), the difference v — u will saturate at t ~ rg, when its magnitude 
is 0{us); for t > Ts, v — u will continue to fluctuate with amplitude us and correlation time ts- 
In this end phase, the equation for Ax will have the form: 

Ax(i) u-(i|x(rs),rs) - u-(i|x-(rs),Ts) (11) 

i.e. Richardson law for the separation of two fluid particles, which, at time rs, were distant 
Ax{ts) = |x(t5) — x''(ts)| ~ S. Of course, when Ax{t) > L (and t > tl), the fluctuations at x 
and xp will decorrelate and Richardson diffusion will turn into normal diffusion. In the extreme 
case of very fast turbulent fluctuations, i.e. tl < ts, the Richardson diffusion phase may be absent 
and separation will cross over directly to a regime of normal diffusion, without passing through 
the regime described by Eqn. (11), with |u — v| saturating to ul instead of to us- 

The picture is complicated a bit by the presence of external forcing on the particle; this leads 
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to new characteristic scales for velocity, length and time: 

UG=9Ts, G = L{gTs/uLf and tq = TL{gTs/uL)^ , (12) 

associated respectively with the drift of the particle in the fluid, and with the size and lifetime 
of the eddies having this characteristic velocity. The relative importance of external forces and of 

random acceleration by turbulence is fixed by the ratio uq/us- Both limit regimes ts — 0,ug ^ 
and Ts > 0, g = are physically realizable, an example of the first one being that of charged 
buoyant particles {p = Po) in a strong electric field. External forcing will have a negligible effect if 
ug < min(us,uz,), i.e., from Eqn. (12), if: 

> —y, — ^^^^ 

mm(l,rs/ri)2 

Thus, inertia may be dominant, when turbulent fluctuations are so strong and fast that ul > grs 
and tl < ts (for p/po ^ 1, in which case external forces are an important factor, this happens 
typically in the vicinity of obstacles) or otherwise, when the particles are small and the Stokes 
time is short enough: ul > g(TsTL)^- 

In the presence of external forces, the difference u — v will fluctuate with amplitude us and 
time scale ts, around uq instead of around zero. For uq ^ us, trajectory separation develops in 
the same sequence of phases as in the case without forces; the time tq < ts identifies only the 
value of the crossover time tb corresponding to |uo — vo| ~ wg- 

When external forcing dominates over the effect of inertia, i.e. when ug ^ us, the sequence 
of phases is altered. The initial phase t < tb remains ballistic in nature. Also the next one, 
for Ts <C t <C rs, is not modified; since gt <C UAa;(t)) the external force contribution is small 
and the particle separation follows the same pattern a la Richardson as in the 5 = case. For 
Ts <^ t <^ TG, however, a new phase will occur: the particle is drifting with velocity ug and the 
contribution to separation from vortices of size Ax{t) is negligible: Ax{t) <C G and u^rc^t) ^ "G- 
The end phase of pure Richardson diffusion described by Eqn. (11) starts only when t > tq, so 
that UAx{t) > Ug, but may be absent in the case of external forces dominating over turbulence: 
Ug > Ul- In this case, also the end phase will be dominated by the drift, with superimposed a 
normal diffusion correction produced by eddies at scales L. While in the absence of gravity, the 
cross-over to fluid-like behaviors occurred at the space scale S, in the present case the transition 
occurs at the scale G, which, since ug ^ us, is larger than S. 

In the three cases of dominant inertia, intermediate (in which external forces dominate inertia. 



but are dominated by turbulence) , and external forces dominating over both inertia and turbulence, 
the separation process will develop in the following possible sequences of phases: 



Richardson {ts < t < tl) Normal {t > t^); 

Intermediate Tg < tq < t^: Ballistic {t < tq) Modified Richardson {tb < t < Tg) Drift 

{ts < t < To) — » Richardson (tq < t < tl) — > Normal {t > r^). 

Dominant External Forces tq > T]^: Ballistic {t < tq) — > Modified Richardson [tq < t < 
Ts) Drift {t>Ts) 

The considerations that have been carried on so far correspond to the intuitive fact that the 
turbulent small scales, i.e. those for which / < niax(5, G), are unable to transport the solid 
particles. For I < S this occurs because of inertia; for < Z < G, it is the external forcing produced 
trajectory crossing which is responsible for this. In the case of the large scales I > max(S', G), the 
contrary occurs: one has |u — v| ^ max(u5,UG') < and the deviation of the particle trajectory 
from that of a fluid parcel on the scale of a vortex of size I will be ^ us.g{ts,g'''i)^ < I and the 
solid particle will behave like a fluid one. 

This scale separation of the turbulent fluctuations and the the subdivision in phases of the 
trajectory separation process, has a counterpart in the self-diffusion (|u(<|xo,0) — uqP), which is 
the relevant quantity to determine the transport properties of the particles [u(t|xo,0) is the fluid 
velocity sampled by the solid particle, which at time t = was in Xq ]. Following we can split 
the velocity difference u(t|xo, 0) — Uq into the different contributions along, and transversal to fluid 
trajectories. The first contribution is just the Lagrangian change u'^(t|xo,0) — Uq; the second one 
u(t|xo,0) — u''(t|xo,0) u^{x{t),t) — u^{x^{t),t) is a complex quantity feeling both effects of 
trajectory separation and of the relative diffusion of fluid elements. In the ballistic regime, using 
Eqn. (3), we have: 



Dominant inertia tq < Tg: Ballistic {t < tb) 



Modifled Richardson {tb < t < ts) 



u(i|xo,0)-u^(i|xo,0) 



u''(xo + voi, t) - u^{ko + uo<, t) 



u 



Uo-Vo 



t ^ ul{\uo - wo\t/L)3 



while, from Eqn. (4): 




Therefore: 



(|u(i|xo,0)-uo|') -GK(e|vo-uo|i)* ^ CK{luAt)i 



(14) 
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where t <€.tb and A = S,G. 

In the opposite regime ts^g^ from Eqn. (11), the transverse separation will be dominated by 
the relative diffusion of fluid elements: u(t|xo, 0) — u''(t|xo, 0) ~ u''(t|x(Ts), rs), — u''(f|x''(Ts), ts) 
while the "longitudinal" component remains u''(t|xo, 0) — Uq ~ (Cpe^)'- Also the first component 
obeys normal scaling, but its presence comes simply from having taken a wrong initial condition for 
u(i|xo,0) ~ u''(t|x(T5), Ts). Said in a different way, the transverse component is not independent 
of the longitudinal one, and they sum together in such a way that the total amplitude is equal to 
that of the longitudinal piece. The total velocity variation is therefore: 

(|u(t|xo,0)-Uo|2) -Coet TL:^f>Ts,TG (15) 

i.e., for long times and turbulence strong enough to prevent drift dominance at large t, the self- 
diffusion properties of fluid and solid particles will be the same. 

In the middle range tb <C f -C ts, all the contributions to u(t|xo, 0) — uq are of the same order; 
these are the Lagrangian change u''(t|xo,0) — uo, the one from Richardson diffusion starting at 
times tb <^t: u^(t|x(rB), tb) — u^(<|x''(tb), tb), and the remnant u(t|x(TB), tb) — u''(f|.-r(rB), tb), 
which is due to trajectory separation. In fact, t identifies a scale I, through the relation tj = i, which 
is the size of the vortices contributing most both to u(<:|xo,0) — Uq, to Richardson diffusion and 
to the trajectory separation |u(f) — v(t)| ^ ui. Even after elimination of the Richardson diffusion 
by choosing a different Lagrangian path, the difference u(f|xo,0) — uq will receive a contribution 
proportional to ui due to the 0{l) displacement between the solid and fluid particle trajectory. This 
is of the same order as the Lagrangian piece u''(f|xo, 0) — uq, associated with the time evolution of 
a size / vortex, which is 0{ui) in the time t ^ Ti. The end result will be again, the normal behavior 
for velocity self-diffusion, described in Eqn. (15), but with a renormalized Co to account for the 
trajectory separation \u{t) — v(<)| ~ ui. 

III. Transport by a synthetic turbulent field 

One of the difHculties in studying the divergence of solid and fluid particles is that the pro- 
cess is affected simultaneously by Eulerian and Lagrangian properties of the turbulent flow. To 
simulate numerically the trajectory of a single solid particle, one would expect therefore that a 
whole advecting fleld should be generated. Clearly, this is as difficult as generating a high Reynolds 
number turbulent flow by means of a DNS (direct numerical simulation) of the Navier-Stokes equa- 
tion. The strategy to bypass this problem, while working in high Reynolds number conditions. 
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has been to generate numerically an advecting flow only in the neighborhood of the solid parti- 
cle. This flow satisfies a restricted set of properties of real turbulence, namely, incompressibility, 
Kolniogorov scaling for the space structure function and normal scaling for the Lagrangian time 
structure function, as described by Eqn. (4). 

The flow was generated by means of a superposition of vortices with the following properties 
(more details are given in the Appendix): 

i The vortices are distributed in a sequence of shells: those in shell n have size / € 
p-n-i 2^"] and are distributed uniformly in the logarithm of scale In this way, 
if + 1 is the total number of shells: < rt < A^, and thus a Kolmogorov length 
rj — 2^^^^ and an effective Reynolds number Re = L/rj ~ 2^+^ are introduced in a 
natural way. 

ii The vortices have a finite lifetime t; — (31^ , with (3 a fixed constant, and are generated 
in the neighborhood of the particle in such a way that the particle lies on the average 
in the support of only one vortex per shell. The radius of this neighborhood scales 
linearly with the shell scale, so that all shells have, on the average, the same number 
of vortices. 

iii A vortex in shell n translates rigidly with speed equal to the value at its center of the 
velocity field produced by all vortices in shells n' < n. 

For the kind of analysis of scaling and crossover behaviors we are interested in, two dimensions are 
expected to be sufficient. In addition, working in 2D allows to neglect difficulties associated with 
the choice of eddy shape, which here is taken as circular. 

The idea of mimicking a turbulent flow with a superposition of rigid vortices, whose dynamics 
imitate as much as possible that of the real ones is not new (see for an application including 
the effect of Lagrangian sweep). In our case, thanks to the limited number of vortices per shell, this 
allows a drastic reduction in the number of degrees of freedom in the system. One more advantage 
is associated with the fact that a flow space structure is present and incompressibility can be 
enforced. This prevents the difficulties present in more standard Lagrangian stochastic models, 
which are associated with the need of a well mixed condition ||ll|] imposed from the outside. The 
same problem occurs in the random walk model considered in in which the vortices are 

modelled as scattering centers for the particles. 



9 



The same approach based on the use of a superposition of vortices localized close to a single 
particle can be utilized to analyze the separation of particle pairs. The separation of trajectories 
in the two cases of a pair of fluid particles (rs = 0) and of a ts/tl = 0.26 solid particle and a fluid 
one are shown in Fig. 1 below. 




Figure 1: Two pairs of trajectories separating in a turbulent flow, in the absence of external forces. 
Effective Reynolds number L/rj ~ lO'^ (17 shells). Trajectories a, b and c refer to fluid particles; 
d is the trajectory of a solid particle with ts/tl = 0.26. Notice the smoother path followed in d. 
Notice also the long transient in the first pair, where both particle are fluid and all the separation 
is due to fluid trajectory divergence. 



In the first case, an initial separation 0.25r/ was adopted; notice, however, the long tran- 
sient associated with the particles lying at "viscous range" separations. The corresponding rel- 
ative diffusion in velocity: (|Aui(i)p) = ([ui(t|xo,0) — wf (i|xo, 0)]^} and in space (|Axi(t)p) = 
([a;i(i|xo, 0) — (t|xo, 0)]^) are shown in Figs. 2 and 3. 

3 

The parameter [3 has been chosen to maximize the ratio C^/Cq, and therefore the velocity 
variation due to trajectory separation. This choice of parameters resulted in: C^^Cq — 2.828. In 
the artificial fiow considered, the turbulent dissipation e is clearly undefined, and the constants Ck 
and Co cannot be determined individually. Assuming Ck = 2, as in three dimensional turbulence, 
would give Cq — 8, which is rather large, as are small the values of the constants c and c entering the 
expressions for velocity and coordinate relative diffusion (|Aui(i)p) — clt and {\Axi{t)\'^) — cet^ . 
One finds: c/Cq — 0.29 and c/Cq, which would lead to: c ~ 2.3 and c ~ 0.025. All of these 
constants are characterized by a strong scatter of experimental values; if the accepted value for Co 

while recent measurements gave c ~ 0.25 



is around 4 |21 , estimates for c went as low as 0.05 
p3[. The uncertainty in c is magnified in c due to the higher power of t, and a factor 10 in c 



10 





11 



corresponds to a factor ^ 2 in c. Nonetheless, the small values of the ratios c/Cq and c/Cq (which, 
contrary to c and c, are independent of Ck) indicate that, although the synthetic turbulence 
field satisfies the proper Lagrangian and Eulerian scaling, the ratio of relative diffusion to self- 
diffusion may be smaller than in real flows. We must not forget, however, that the relevance to 

3 

three dimensional turbulence of parameters such as e.g. C^/Cq, obtained from a two dimensional 
synthetic field, may be somewhat limited. 

The analysis of trajectory separation, contrary to that of the self-diffusion (|u(t|xo,0) — uqP) 
is plagued by several difficulties. The most important one is the limitation in the sample size: 
small and large separations, in the relative diffusion case, are calculated from the same sample 
of trajectory pairs, while, in the self-diffusion case, the sample size is inversely proportional to 
the scale in exam. This limitation is particularly serious, since most of the information on the 
separation of solid particle trajectories is conditioned to the initial velocity difference. The second 
difficulty is the strong effect of transients occurring when the particles are very close and stick 
together, which is absent in the case of self-diffusion. In fact, as suggested in |Q, to obtain the 
curves in Figs. 2 and 3, it was necessary to average over the logarithms of the doubling times of 
|u(t|xo,0) — u''(t|xo, 0)1 and |x(i|xo,0) — x''(i|xo, 0)|, instead of at fixed times, over |u(i|xo,0) — 
u^(i|xo,0)| and |x(t|xo, 0) - x^(<|xo, 0)|. 

The self-diffusion properties of particles with different values of the ratio ts /tl are shown in 
Figs. 4, and 5 for the two cases without and with external forces. 

The time-scale of fluctuation for the difference u — v appears to be of the order of tq and ts , 
respectively, in the forced and unforced case. In the unforced case, the largest ts/tl — 2.6 curve 
is basically indistinguishable from what would result from setting v = 0, i.e. ts = oo. (The short 
time behavior of this curve, incidently, gives then also an idea of the applicability of the Taylor 
frozen fleld hypothesis in the absence of a mean flow). In the forced case, only the largest value of 
the drift ug leads to an appreciable decrease in the correlation time for the fluid velocity sampled 
by the particle. In the other cases, in which > uq, the synthetic turbulence model is probably 
too crude to give any definite answer. 

The inertial range is compressed passing from space to time (compare with Fig. 11 in the 
Appendix) and this has the consequence that crossover behaviors dominate also at very high 
Reynolds number. In order to observe cleaner scaling behaviors in the velocity self-diffusion, it 
would be necessary to condition the averages to the initial value of the difference u — v, as shown 
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10- 

10-^ 10-^ 10-2 10-1 IQO IQi 

Figure 4: Velocity self-difFusion in the inertia dominated ease {ug = 0) for different values of 
ts/tl- ' — I — ts/tl = 2.6; heavy line: ts/tl = 0.26; ' — o — ': ts/tl = 0.026; thin line: ts/tl = 
0.0026; dots: fluid particle. The Kolmogorov time-scale is t^/tl ~ 3.8 x 10"^, corresponding 
to an effective Reynolds number Re ~ 2 x 10^ (18 shells). In the various cases, we have, for 
<-v| = vp) -^i^: ts/tl = 0.0026: mZ'<_v| = 0-0049; ts/tl = 0.026: u^V^,,^, = 0.047; 
ts/tl = 0.26: wZ'<-v| = 0-39; ts/tl = 2.6: ^Z'^-vl = 0-9. 




Figure 5: Velocity self-diffusion in the external forcing dominated regime (rs = 0) for different 

values of ug/ul] same flow parameters as in Fig. 4. ' — \ — ': ug/ul = 1.83: heavy line: ug/ul = 
0.09; ' — o — ': ug/ul = 0.01; dots: imforced case. In the various cases we have: uq/ul = 1.83: 
«Z'<-v| = 1-7; ug/ul = 0.09: wZVf„_^| = 0.11; ug/ul = 0.01: mZV^^,^, = 0.015; 
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in Fig. 6 in the unforced case. 




10-^ 



10-^ 



10-3 



10-2 



10-1 



10° 



IQi 



t/TL 



Figure 6: Velocity self-difFusion conditioned to different values of the ratio (|u — v|)-i|u — v|, for 
fixed ts/tl = 0.26, ug = 0. Same flow characteristics as in Fig. 4. ' + ': 2 < (|u — v|)-^|u — v|. 
' □ ': (|u — v|)~^ |u — v| < 0.12. ' O ': unconditioned. Notice how ballistic behaviors are enhanced 
in the two conditioned cases. 

An example of cross-over behavior is the rcnormalization of the constant Co over the whole 
inertial range for values of ts/tl as small as 0.0026. 

It is possible to have an idea of what is happening, at least in the unforced case, focusing on 
the separation of trajectories in a vortex of size I > S. The separation process is the result of the 
particle being accelerated by vortices of all the available scales Prom Eqn. (1), large vortices 
with V > I will produce a constant velocity difference ui'Ts/ti'. The contribution to trajectory 
separation by a single vortex of size I', will be, in the lifetime Tf. 



The contribution from vortices at scale V G [S, I], instead, is similar to a random walk, since each 
vortex I' produces an independent displacement of the particle in the time Tf, the total displacement 
produced in this time will be therefore: 



Smaller vortices, of size I' < S, produce altogether an 0{S) contribution to displacement. Exploit- 
ing the fact that the vortices are distributed uniformly in the logarithm of scale, and combining 



Al{l') ~ TiUpTs/Ti: ~ Its/ti{ui/ui>). 



(16) 



Al{l') ~ 1'ts/ti'{ti/tv)^^ ~ Its In 



(17) 
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with Eqns. (16-17), the velocity shift wiU be: 

S<1'<1 S<1'<1 S<1"<1' ' 

The effect on the time velocity structure function will be, from Eqn. (6): 

(|u(t|xo,0) - uoH ^ Coetil + ^[\og{bt/Ts)r) (19) 
with a and b constants. A fit of (|u(i|xo,0) — uqP) using Eqn. (19) is shown in Fig. 7. It is to 
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Figure 7: Negative renormalization of Co for small values of ts/tl; O: fluid particle, fitted by 
Eqn. (14) with Co = 8 (for Cr = 2). +: ts/tl = 0.0026 ug = fitted by Eqn. (11) with Co = 8, 
a = —0.3 and b = 1.4. Using still Eqn. (14), this would correspond to a renormalization Co — 7. 

be noticed that the renormalization of Co is negative in all the inertia dominated situations, with 
the exception of the T5 > 1 cases, in which the linear spectrum is essentially absent. A negative 
renormalization of Co corresponds to an inhibition of the time variation for u(i), which is akin to 
the increase of the correlation time for u(t), predicted in ^ in the case of a purely large scale 
flow. It is not clear, however, how much these phenomena depend upon modelling assumptions. 

The crucial aspect seems to be whether the time-statistics is given in an Eulerian frame or 
along trajectories. An Eulerian time-statistics would suggest a faster fluid velocity variation in a 
Lagrangian frame p5| , which is the physical motivation of the results obtained in ^ . Following 
the same argument, a time-statistics given along trajectories would lead to the fixed particles 
feeling the fastest velocity variation. Analytical calculations for random velocity fields with short 
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correlation time and Lagrangian time statistics, however, confirm the numerical result of slower 
velocity variations along particle trajectories p^ . 

IV. Application to Lagrangian transport modelling 

The analysis carried on so far indicates that the strongest effect of particle inertia on velocity 
self-correlation is an enhancement of fluctuations at time-scales below Tg and tq. Contributions 
at long time-scales, associated with renormalization of the constant Cq, remain present even for 
very light particles, but the correction is small compared to the uncertainties in the value of this 
constant, even in the case of fluid trajectories. Hence, for time scales above Tg and tq, it is 
acceptable to set u ~ u'' and to model the fluid velocity u sampled by the heavy particle, with the 
same Langevin equation used for u'', i.e. the second of Eqn. (2). Being able to model the difference 
u — u'' becomes necessary, however, when ts^g ~ tl and trajectory properties at time-scales below 
Ts and TG become important. 

In a Markovian approach like that of 0, the equation for u would be obtained from the 
one for u'', by adding, at discrete times tn — riA, A < ts.Gi uncorrelated increments Au extracted 
from a PDF (probability distribution function) in the form P(Au| Ax, u'^), with Ax ^ (v — u)A 
the trajectory separation in the time A. The shift Au will have a mean part Au(u''), which is zero 
when u'' = and a fluctuation Au: Au ~ u/^x The cumulative effect of these increments would 
be felt also at time-scales above tq and r^, in the form of additional noise and drift terms in the 
equation for u. In particular, the amplitude of the noise correction is estimated from |u — v| ~ us,g 
and Eqns. (3,8,12): 

A-i(Au)2^e-(rs,G/A)^ (20) 

which is an unphysical result: this amplitude remains finite also in the fiuid particle limit ts^g I^l —>- 
0. The most the drift term can do to improve upon the situation, is to make sure that (u^) — u\ 
and that the only corrections be in the value of the constant Co and the Lagrangian correlation 
time tl seen by the solid particle. Neglecting the drift, which amounts to extracting the velocity 
increments from a PDF not conditioned to u'', would lead to {v?) > m|; this would be equivalent 
to say that solid particles privilege higher velocity regions of the flow, independently of the value of 
Ts.G- Going back to the discussion preceeding Eqn. (15), neglecting the conditioning on u'' would 
mean that any correlation among velocity increments along fiuid trajectories (longitudinal) and 
from trajectory separation (transverse), is disregarded. (It is to be mentioned that, contrary to 
the prediction of |^], the renormalization of the turbulent correlation time, necessary to the fluid 
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velocity amplitude constant, leads to a decrease rather than to an increment of tl). Notice that 
approximating spatial correlations with exponentials pO| , is not sufficient to cure this problem; 
the only gain is eliminating the spurious noise dependence on the discretization A. In fact, if we 
neglected Kolmogorov scaling of space correlations and approximated Au u\l/A, we would still 
have A^^(Au)^ ^ u\/ta ^ e, A — S,G, which remains finite when ts.g ~^ [compare with Eqn. 
(20)]. In other words, anomalous diffusion is not an essential ingredient; memory effects are. 

Eliminating cumulative effects of the increments in u leads automatically to the introduction of 
memory and to a non-Markovian model. A simple Gaussian model with these properties is obtained 
proceeding in the same philosophy of section III: the velocity u is modelled by a superposition of 
solutions of stochastic linear equations, each associated with a spatial scale of the turbulent ffow. 
For homogeneous isotropic turbulence (summation over repeated vectorial indices is understood): 

< M_^-iA^.(a«)^(£^)^e.(i) (21) 

Mt)i,m - 5.,r-iS(t/f„); /dxS(a;) = 1 
where A" — —Ml^uJ, f„ and M^" are function of u — v, c„ ~ 2 i~Co/2, and for |u — v| = 0: 
f„ = 2^"fo. The noise ^i{t) must have finite correlation time in order to avoid that the dynamics 
of u at a given time-scale r receive contribution by shells with f„ 3> t (this is implemented 
automatically in a numerical code, using timesteps proportional to f„, in the different shells). 
With these prescriptions, for |u — v| = 0, u" will ffuctuate on time scale f„ with amplitude 
(|u"|2) = 2-"(|u"|2) oc f„, and for t < tl Eqn. (15) wiU be satisfied. 

The variable u"(t) in Eqn. (21) represents the contribution to the Lagrangian velocity com- 
ponent u{t) from vortices of size l{n) such that ti(^„-j ~ f„; hence l{n)/L ^ and then also 
(|u"|^) ^ ^'[{n) UL{l{n) / L)^ . If |u — v| ^ u" is the frozen velocity field of the vortex of 
size l{n). We have therefore: 

T^n _ J 1 |U - v| < 



T„ 



(22) 



t(0) [7Ui(„)/|u - v| |u-v|>Ui(„) 
where f„(0) indicates the value of f„ for |u — v| =0 and 7 is an 0(1) constant. For |u — 
v| ^ ui(n), ('""(^)'"j (0)) is just the space correlation function at separation (v — u.)t of the flow 
field produced by vortices at scale l{n). Incompressibility fixes then its tensorial structure Q: 
{uf {t)uj (0)) = f{t/fn)Lij + g{t/fn){5ij — Lij) with Lij — ^" l u'-l^:^ s-nd in three dimensions: 
g{x) = f{x) + ^f'{x). Using Eqn. (13) also for |u— v| ^ is equivalent to assuming exponential 
space correlations, and from f{x) = exp(— x) we find g{x) = (1 — |-)exp(a;) ~ exp(— Ix). In the 
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opposite limit |u — v|/u;(„) — * 0, isotropy imposes (u"(t)u"(0)) = f{t/Tn)5ij. These conditions fix 
the form of the tensor M"j in the two opposite regimes: 



|U - v| < Ui(„) 

^Lij |u - v| > u/(„) 



(23) 



The model described by Eqns. (21-23) is reminiscent of a continuous time random walk |27 



that the non-Markovian nature of the process originates from varying the correlation times of the 
various contribution to u(t), rather than their amplitudes, which remain invariant. 

The self-diffusion properties of the solution to Eqns. (21-23), limited to one dimension, are 
compared in Figs. 8 and 9, with the synthetic turbulence results of section III. The time tl and 
the noise amplitude Cqc in Eqn. (15) have been selected to fit the synthetic turbulence result 
for Ts = tg = and the expression f„/f„(0) = (1 + 0.3|u — v|/u/(„))^-^ has been used for the 
correlation times. When f„ < fN{0), with N the fastest shell, mode n was disregarded and the 
fastest mode with correlation time above 2~^tl was renormalized to account for the energy of the 
lost shell. 



s 



s 
I 




Figure 8: Comparison of solution to 1-dimensional version of Eqns. (21-23) (continuous lines) and 
synthetic turbulence results for velocity self-diffusion. 'O' ts/tl — 2.6. ' + ' ts/tl = 0.26. '□' 
ts/tl — 0.026. ' X ' = 0. Equation (21) has been solved with 13 modes and fo(0) ~ 0.35rL. 



The model just described allows to keep into account the Kolmogorov-like nature of the small 
scale spatial structure of high Reynolds number turbulence. It is easy to see this. From the 
discussion in Section II, the typical particle-fluid relative velocity is |u — v| ^ max(us,UG)- The 
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Figure 9: Another comparison of solutions of Eqns. (21-23) (continuous line) and synthetic turbu- 
lence results. Here, self-diffusion is conditioned to different values of the ratio (|u — v|)-i|u — v|, 
for fixed ts/tl = 0.26, ug = 0: ' + ': 2 < (|u - v|)-i|u - v|. (|u - v|)-i |u - v| < 0.12. 'O': 

unconditioned. Equation (21) has been solved using the same parameters of Fig. 8. 

fluid velocity increment along a solid particle trajectory (|u(i) — u(O)p) will be dominated by that 
shell for which fn{us,G/ui(^n)) — t ^-nd t — ts^g will give the cross-over from the situation in which 
the corresponding f„ obeys the first or the second of Eqn. (22). From Eqn. (22), and the scaling 
for the shell characteristic time and velocity: f„(0) = to2-", (|u"p) ^ uj^2^"~^, we thus obtain: 

At short times, the fluid velocity self-diffusion will have the scaling, which comes from the 
spatial profile of the small vortices being seen as frozen by the crossing solid particle. Conversely, 
at long times, we will find the usual normal scaling of the Lagrangian time structure function 
[see Eqn. (4)], which is obtained disregarding trajectory separation, and equating t with the eddy 
turn-over time. 

The problem of Markovian models |l^ , discussed at the beginning of this section, disappears 
in the present approach. Trajectory crossing effects, in the limit ts^g ^ 0, are confined to shells 
with n — > oo, and their contribution to T]^ disappears in the limit (as it should). In this respect, 
this model is closer to the ones by ^ , in which the effect of particle inertia is accounted for by 
correlation time rcnormalization, rather than modelling directly the change of fluid velocity from 
trajectory crossing. Of course, the intent of modelling the turbulent small scales, adopting either 
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the Markovian approach of |l^, or the non-Markovian approach used here and in jl^. 15 , was 
completely absent in p8{ . 

In basically all existing theories of heavy particle dispersion, the Kolmogorov structure of small 
scale fluctuations is disregarded. Typically, an exponential |l^, |l^, ^ or Gaussian |j profile for 
the space correlations is adopted. In some approaches ^, H, also the Lagrangian nature at short 
time separations, of the time correlations, is neglected. 

After first submission of the present work, this author became aware of the paper by Reynolds 
and Cohen p9[ | , in which a fractional Brownian motion algorithm was used to generate the anoma- 
lous scaling behaviors resulting from trajectory crossing. The advantage of the present approach 
is the predictive ability for the scaling, whose physical mechanis is built in the algorithm by means 
of the shell superposition technique. It is interesting to analyze the prediction of the present 
model on the test case considered in based on the Snyder and Lumley experimental data . 
Figure 10 illustrates a power law fit for the fluid velocity self diffusion along a solid particle trajec- 
tory simulated with Eqns. (21-23), with particle parameters corresponding to the "glass", "corn 
pollen" and "hollow glass" particle case considered in pO]. The same form for fVi used before: 



S 
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Figure 10: Fluid velocity self-diffusion predicted by Eqns. (21-23) in three test cases; □ (glass): 
ug/ul — 1.76, ts/tl — 0.26; O (corn pollen) ug/ul — 0.7, ts/tl — 0.1; -i- (hollow glass) 
ug/ul — 0.052, ts/tl — 0.01. The continuous line illustrate the power law fit by t^^ \ in the three 
cases: H ~ 0.41; H = 0.45 and H ~ 0.5. 



T„/r„(0) = (1 + 0.3|u — v|/u;(„-)) ^, leading to the results of Figs. 8 and 9, continued to be used. 
In the three cases the fit gives for the exponent H: (|u(i|xo, 0) — uqP) ^ t^^ , H ~ 0.41, H ~ 0.45 
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and H ~ 0.5. This agrees well with the results in with the exception of the hollow glass case, 
in which a behavior closer to that of a fluid parcel is predicted here. 

A shortcome of the present model is clearly that the exact form of the shell time-scales must be 
provided from the outside, since Eqn. (22) is only valid asymptotically. Among the other things, 
knowing the exact form of T„/f„(0) for |u — v|/u;(„) < 1 would be necessary, at least for small n, to 
determine the fluid velocity correlation sampled by the solid particle. This shortcome is common 
also to other models ^ , in which the renormalization of the correlation time by inertia must be 
introduced as an external ingredient. It is to be stressed, however, that the inertia renormalization 
of the correlation time belongs, together with the determination of the constant Co and of the 
relative size of Eulerian and Lagrangian correlation times , to a class of problems which are 
havily influenced by the turbulent large scales. It is thus to be expected that this parameter be 
dependent on the nature of the flow and it be necessary to provide it from the outside. 

It is possible to extend the Lagrangian modelling approach illustrated in this paper to non- 
Gaussian inhomogeneous turbulence situations, imposing, in the fluid particle limit |u — v| = 0, 
individually on each shell, the well-mixedness condition It must be stressed that the well 

mixedness condition can be imposed only in the fluid particle limit v = u. Apart of all the 
problems connected with ergodicity of solid particle trajectories, the particle concentration profile 
in the inhomogeneous turbulence flow is unknown and, in fact, is typically one of the desired 
outputs of the model. 

Basically, one has a condition on e(x, t) from the integral time scale: 

Gotl(x, t) 

plus a condition for each of the (tensor) moments which are experimentally available: 

5]([uT)(x,i) = (u^')(x,i) 

n 

with Kolmogorov scaling giving, for p — 2: (u^)(x, i) = 2^"^^m|(x, t). The choice of the contri- 
bution from each shell to the non-Gaussian part of the statistics remains arbitrary, conditioned to 
statistical realizability of the moments for the Un- 

A possibility could be to partition the contribution equally among different shells: ([u"]^')c(x, t) = 
2~^^ ' {[u]P)c{x,t), with subscript c staying for cumulant. A choice which is more consistent 
with the standard practice in stochastic Lagrangian models, to let the non-Gaussianity be gen- 
erated by a non-linearity in the drift term, considering it therefore to be a purely large scale 
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feature of the flow, would be to concentrate all the non- Gaussian! ty production in the first shell: 
(|u"|P),(x,t) = {uP),{^,t)di,„. 

Once the PDF's P„(u",x, for the various shells are obtained from the cumulants through 
some reconstruction technique the Fokker-Planck equation for the second of Eqn. (21) can 
be inverted to determine the drift terms A" from e and P„; considering for simplicity the case of 
stationary turbulence: 

As it is well known (see e.g. |3^), in inhomogeneous turbulence situations and for (ti") — 0, Eqn. 
(25) leads to a constant contribution to the drift, which counterbalances the drift A" from the 
gradient term u'ldiPn and enforces zero mean for u". 

Notice that this balance disappears in the case of the solid particle flow, in which case the 
particle velocity will be characterized by a non-zero turboforetic mean drift. The multivariate 
PDF P({u"}, v,x) corresponding to a uniform particle concentration will obey the Fokker-Planck 
equation: 

vAP - a.,( ^' ~"'~"'^' P) + V d^Af-'A-P) = V £!l!l^i9^„a,.(P/f„) (26) 

TS — ' — I ' ' 

n n 

from which, for small ts/tl (and then us <^ u^), the drift can easily be calculated. Multiplying 
Eqn. (26) by Vi, integrating over v and {u"}, and using {viVj) = (uiUj) + 0(u|) — (uiUj), we 
obtain in fact: 

{vi) ^ -Tsdj{uiU;j) (27) 

Thanks to the well-mixedness condition, however, this turboforetic drift is free of spurious contri- 
butions from unproper Lagrangian modelling of the fluid limit case. 

It must be said that, contrary to what happens (at least in part) in a model like Jill , the time 
structure of the individual trajectories is here physically relevant. This because of the coupling 
between u and v through Eqn. (1). Trying to base an extension of the present model to non- 
homogeneous non-Gaussian turbulence, on a statistical constraint like well-mixedness, may thus be 
of limited value. It must also be said that the ergodic property is not satisfied in general by heavy 
particles. (In the case of fluid elements, the ergodic property is always assumed, implicitly, as a 
working hypothesis). It is not clear, therefore, whether it is appropriate to hypothesize invariance 
even only of the amplitude for u(<), when passing from |u — v| = to the heavy particle case 
|u- v| > 
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V. Conclusions 

The analysis carried out in this paper provides a picture of the behavior of heavy particles 
in turbulent flows, and, in particular, a prediction for the scaling properties of the fluid velocity 
sampled by the particles, in function of the flow and particle properties. The predictions based on 
scaling analysis were confirmed by synthetic turbulence simulations and allowed the derivation of 
a Lagrangian transport model taking into account the non-Markovian nature of the fluid velocity 
along the heavy particle trajectory. 

It is important to understand that this non-Markovian nature cannot be overlooked without 
getting into contradictions. For instance, it would be rather natural, in a model like to 
hypothesize a normal diffusion behavior for u'' at high Reynolds numbers, and to study pairs 
of particle-fluid trajectories of duration A < rs, in which the approximation u ~ u'' is used. 
But, as discussed at the beginning of Section IV, if one neglects correlations between the velocity 
increments Au = u — u'' at the discrete times i„ = nA, modiflcations in the time correlations 
sampled by the solid particle will appear up to time-scale tl and remain present in the fluid 
particle limit ts.g ^ 0, which is unphysical. Even worse, if the PDF P(Au|u''), from which Au is 
extracted, is not chosen carefully enough, and the condition on is disregarded by choosing say 
P(Au) = P(u'=(x + Ax, t) - u'=(x, t)), one will flnd also (uf ) 7^ ul- 

The scaling analysis suggests a picture of transients between different spectral slopes at different 
time-scales, in which the controlling parameters are rs and us {tq and uq in the case of strong 
external forces), which are the time and velocity scales for u — v. Both synthetic turbulence 
simulations and the runs from the Lagrangian model, however, indicated that cross-over behaviors 
rather than pure scaling dominate even for very high Reynolds numbers. 

In the last section, some ideas for inclusion of non-Markovian effects in a Lagrangian stochastic 
model for solid particle transport have been illustrated. In this way, the Kolmogorov nature of the 
turbulent small scales and the Lagrangian character of time correlations is fully taken into account. 
The suggested approach is based on the use of a superposition of Markovian processes, which will 
be slower than a simpler Markovian model by a factor of the order of the number of terms in the 
superposition. (In our case, with 12 shells, the ratio was approximately 8). A similar approach 
has been recently used to model turbulence intermittency 

An alternative approach to generate the sub-diffusive behavior for the velocity, at time-scales 
below ts^Gj could have been to resort to a continuous time random walk |^, or, equivalently. 
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to Levy walks p3[ . This would have been lighter from the computational point of view, but the 
control on the form of the PDF for the velocity and the velocity difference would have been much 
more difhcult. The control on the shape of the correlations for long times, would have been difficult 
as well. Notice that, in the present approach, recirculation phenomena and correlations modelled 
on Frenkiel functions ^ could easily have been taken into account in Eqn. (21) by the use of 
complex relaxation times. 

Compared to the recent model by Reynolds and Cohen | |29| , the present approach provides in 
self-consistent way the scaling exponent for the fluid velocity self-diffusion. This both in the case of 
dominant inertia and dominant external forces and with no special cost in terms of simplicity and of 
the possibility of imposing a well mixedness condition in inhomogeneous non-Gaussian turbulence 
situations. 

To obtain a Lagrangian model of practical use, a precise form for the dependence of the shell 
correlation time f„ on u — v should be provided. The two-dimensional synthetic turbulence simu- 
lations, which have been carried out, can at most confirm scaling behaviors and perhaps give the 
sign of the inertia produced correction to the correlation time [see Eqn. (19)]. The form of the 
correlation times f„ could probably be obtained only by DNS of the full Navier-Stokes equation in 
the proper tree-dimensional flow conditions. 

However, although other such attempts to take into account the Kolmogorov structure of the 
turbulent small scales, were already carried on |l^, it is still not clear how relevant these aspects 
are to the derivation of a dispersion model of practical use. Probably, the main interest remains 
conceptual. Considering the reduced information on the turbulent flows available in practical 
situations, and the fact that dispersion is governed by the local large scales anyway, a minimal 
implementation of the model could disregard the small scales altogether and it could be sufflcient 
to keep only one shell in Eqn. (21). In such a hyper-simplified version, the present model would 



become equivalent to 1 28 , with a generalized form of the Lagrangian correlation profile for ts — 



to allow enforcement of the well mixedness condition. 
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A. Appendix 

In this appendix, a closer look is given to the synthetic turbulence model described in section 



24 



III. In this approach, the effect of a homogeneous turbulent flow is modelled by a velocity field 
localized around the particles that one is tracking. This limits the evaluation of spatial objects like 
e.g. correlation functions, to the case in which the correlation points coincide with the particles. 
This approach, has been used to analyze the process of relative diffusion (sec Figs. 2 and 3), but, 
as already discussed, is plagued by difRculties. The space structure functions can be obtained in 
a simpler way from the assumption that the local field around the particle is representative of the 
whole flow. The Eulerian velocity field is in the form: 

u^(x,f) = ^u,(x,i) (Al) 

i 

where i labels the vortex. The contribution from vortex i is 

u,(x,i) = ar/ W^,R(— ^)T(-^) (A2) 

where W, are uncorrected zero mean and unitary variance random variables, a is a constant and 
5ti{t) and ti indicate respectively the vortex center position at time t and the time of "birth" of 
vortex i. We have considered vortices with compact support in space and time; precisely: 

fR(x) = 83 X xi?(l - x) 

\T{t)^H{t{l-t)) ^ ' 

where H{x) = x6(x) with 0{x) the Heaviside step function. The lifetime of the vortex has been 
taken as a deterministic function of the radius: = /3r3. 

Since the W, are uncorrelated, the two-point correlation reduces to a sum of contributions from 
the individual vortices: 

(u(x,t)u(0,0)) = / dlogr ( di ( d2a;iV(r)r3T(^)r(^)R(^i^)R(--) (A4) 

J J J Tr Tr r r 

If any point feels on the average the effect of a fixed number of vortices per shell, the vortex density 

N{r) can be written in the form 

iV(r) =iV(7rr2f^)-i (A5) 

with N a constant, which is exactly l/log2 in our case of one vortex per octave contributing to 
velocity at a given time and position. 

Substituting Eqns. (A3) and (A5) into Eqn. (A5), one can find the expression for the longitu- 
dinal strucure function S2{1) = ((ui(Zei,0) — Mi(0, 0))^): 

^'^^^ = 15^ / 5 / d's^i(x)[i?i(x + (Vr-)ei) - i?i(x)] (A6) 
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Integrating numerically, one finds: 82(1) — O.OOSQa^Zs , i.e. in terms of Kolmogorov constant and 
mean dissipation: 

a ~ 16C|-e3 (A7) 
In Fig. 11 below is illustrated the structure function for an 18-shell synthetic tubulence field. All 
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Figure 11: Space structure function from an 18-shell synthetic turbulence field with a = 22.627, 
L = 1 and fit by a Kolmogorov spectrum S2{x) = 2xs . 

simulations have been carried out with (3 ~ 8.3/a, corresponding to an eddy life of the same order 
of the eddy turn-over time. (Large and small values of /3 correspond respectively to long and short 
vortex lives). 

Once we know the coefficients for x~iS2{x) and t~^{\ui{t) — ■ui(O)p), from Eqns. (10) and 

_ 3 

(AT), we can eliminate the dissipation e to find the ratio Cq ^Ck- We obtain from the simulations 

_ 3 

Cj^^Cq ~ 2.828. If we take for the Kolmogorov constant Ck — 2, we will have Co — 8, which 
is somewhat large. The explanation of this fact may lie in the way a particle deviates from the 
flow line of a given vortex due to the effect of vortices at other scales in the same region. In 
such a model, this effect is magnified by the rigid structure of the vortices, and a particle will 
move out of the one in which it lies at a given time, faster than it would do if the vortex were 
able to be deformed by the flow. More important, however, is probably the fact the this flow is 
two-dimensional and it would be interesting therefore to compare it with a closer physical situation 
like two-dimensional turbulence in the inverse cascade range. 

Following the same approach, it is possible to calculate ul and, from Eqns. (3-4) also L and 

26 



tl- The main parameters of the synthetic turbulent field are listed in Table I below. 



UL ~ 0.61 
L -0.17 

TL ^ 0.2 



Tl ~ 0.07 

C'^Co 2.83 
Cq^c ~ 0.29 



Cq-'c ~ 0.003 

a = 22.627 
P = 0.367 



Tabic I 

The last issue is how the vortices were generated in the synthetic turbulence field. The total 
number of vortices in shell n was taken equal to the integer part of a variable TV, obeying an 
equation in the form: 

J^n = -T-\\u - v|)sign(7V„ - 1) (A8) 

where 7V„ is the number of vortices having the particle in their support at the given time and 
T„(|u — v|) is the average time a particle travelling across the fluid with relative velocity v — u 
feels the effect of a vortex in shell n: 

r„(|u-v|)~2-*F(2-3|u-v|) (A9) 

with F{x) oc for x large. In this way, the number of vortices could be adjusted, in a time of 
the order of the vortex life, in such a way that the particle felt the effect of only one vortex per 
shell. 

If the number of vortices alive at a given time was less than A/", a new one was put in the shell 
in an appropriate domain around the particle and with an age chosen at random between zero and 
the vortex lifetime. The domain was an annular region comprised between r and jr with 7 > 1, 
in such a way that the velocity at the particle position did not suffer from discontinuities. In all 
simulations, we set 7 = 2.3; with this choice, it was enough to keep A/' < 4 in all simulations, 
with the exception of those with ts/tl = 2.6 which required A/" < 5. As in other shell models, the 
main limitation was not the number of degrees of freedom, but the strongly different time scales 
involved when the number of shells becomes large. 

Apart of the potential increase in the number of eddies required in the neighborhood of a 
particle, and of the conceptual problems associated with the eddy geometry, extension of the 
present technique to three dimensions is clearly straightforward. 
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